Pressure-Induced Anomalous Phase Transitions and Colossal Enhancement of 

Piezoelectricity in PbTi0 3 
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We find an unexpected tetragonal-to-monoclinic-to-rhombohedral-to-cubic phase transition se- 
quence induced by pressure, and a morphotropic phase boundary in a pure compound using first- 
principles calculations. Huge dielectric and piezoelectric coupling constants occur in the transition 
regions, comparable to those observed in the new complex single-crystal solid-solution piezoelectrics 
such as Pb(Mg 1/ / 3 Nb 2 /3)03-PbTi03, which are expected to revolutionize electromechanical applica- 
tions. Our results show that morphotropic phase boundaries and giant piezoelectric effects do not 
require intrinsic disorder, and open the possibility of studying this effect in simple systems. 

PACS numbers: 77.80.-e, 77.84.-s, 63.20.-e 
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The classic ferroelectric PbTi03 has been known to 
have a single ferroelectric tetragonal (T) to paraelectric 
cubic phase transition with increased temperature Q or 
pressure Q since its discovery Q in 1950. It seemed 
unlikely to discover any new ferroelectric transitions in 
PbTiOa, but we predict new ferroelectric phases under 
pressure. Piezoelectrics convert electric energy to me- 
chanical energy, and vice versa. They are widely used in 
medical imaging, acoustic sensors and transducers, actu- 
ators, eic Jjjjd- PbTiC>3 has been extensively studied 
QllHllllllEHi3 to understand the electronic ori- 
gin of ferroelectricity. The Pb-0 and Ti-0 hybridization 
in PbTiOa weakens the short-range repulsions and gives 
rise to the ferroelectric distortion 0- It has been found 
that pressure suppresses ferroelectricity since compres- 
sion favours short-range repulsion, but there has been no 
previous work on the effect of pressure on piezoelectricity. 

PbTiC>3 has a high c-axis strain of 6.5% (7.1%) at room 
(low) temperature. Under ambient pressure it undergoes 
a first-order phase transition at T c — 763 K, and this 
transition is regarded as a typical displacive transition, 
which is associated with soft-modes [l|. The displacive 
phase transition temperature T c decreases under hydro- 
static pressure. A Raman study of PbTiC>3 shows that 
T c reduces to 300 K and the phase transition becomes 
second-order at P — 12.1 GPa 0]. It indicates a tri- 
critical point in the phase diagram, where the hrst-order 
phase transition changes to second-order, and it has been 
identified at P = 1.75 GPa, T = 649 K Q. Increasing 
pressure further reduces T c . 

PbTi0 3 is an end member of Pb(Mg 1/3 Nb 2 /3)03- 
PbTiOa (PMN-PT) and Pb(Zn 1/3 Nb2/3)0 3 -PbTi03 
(PZN-PT), which have piezoelectric coefficients an order 
of magnitude larger than those of conventional ferroelec- 
tric simple compounds |l3|. It is also an end member 
of PbZr0 3 -PbTi0 3 (PZT), the most widely used cur- 
rent piezoelectric, which is ubiquitous in modern tech- 
nology. These materials have a common feature, namely 
the morphotropic phase boundary (MPB), and they have 
optimum piezoelectric efficiency near the MPB. Polar- 



ization rotation [tO, Il5| is believed to play an impor- 
tant role for this extraordinary property. Close to the 
MPB, the energy surface for polarization is very flat so 
that for example a polarization along [111] direction in 
the rhombohedral (R) phase can be easily rotated to- 
ward the tetragonal polarization direction by applying 
an electric field along [001] direction, and the R phase is 
transformed to the T phase via the intermediate phase(s). 
This mechanism has been elucidated by theoretical cal- 
culations on BaTiOa and PZT [lHE3, and by ex- 
perimental findings of low symmetry phases [ljj LUi uil • 
The giant piezoelectric effects always do not occur along 
the spontaneous polarization direction (ljj. Il5j|). At am- 
bient pressure PbTiC>3 has no such MPB because it has 
a rather stiff energy surface near the T phase and it has 
no R phase. Under high pressures the energy (enthalpy) 
surface is expected to be softer, and an MPB could arise. 

In this Letter we address both issues of phase transi- 
tions and piezoelectricity of PbTiOa under hydrostatic 
pressures. We perform total energy as well as linear 
response computations, which have been proved to be 
highly reliable for ground state properties. We find 
anomalous phase transitions and giant enhancement of 
dielectricity and piezoelectricity near the phase transi- 
tion regions induced by pressure. 

We have carried out zero temperature ab initio compu- 
tations based on density functional theory (DFT) within 
the local density approximation (LDA). We used the 
pseudopotential planewave method implemented in the 
ABINIT package [2lJ ■ The planewave energy cutoff is 60 
Hartree, and the k-point mesh for Brillouin Zone integra- 
tion is of 6 x 6 x 6. We used the OPIUM program 22] to 
generate norm-conserving pseudopotentials, which were 
rigorously tested against the full-potential linearized aug- 
mented planewave (LAPW) method 23]. We included 
semi-core states of Pb 5d 10 , Ti 3s 2 3p 6 4d 2 , and O 2s 2 in 
valence states. We chose the LDA instead of the gener- 
alized gradient approximation (GGA) because the GGA 
catastrophically overestimates both equilibrium volume 
and strain for tetragonal PbTiOs |24j. On the other 



TABLE I: Elastic (c M „) and piezoelectric (e«/) constants of 
tetragonal PAmm PbTiOs with lattice constants a — 3.902 A, 
c = 4.155 A. Cfj, v and ei„ are in GPa and C/m 2 , respectively. 
Here FS refers to the finite strain method. 
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hand, the LDA moderately underestimates the volume 
(60.38 A 3 ) and strain (4.6%), and the experimental Vq 
corresponds to a negative pressure Pq — —2.2 GPa. 

To calculate the dielectric susceptibility Xy, elastic 
constants , and piezoelectric stress constants e,„ (Here 
Latin indexes run from 1 to 3, and Greek indexes from 1 
to 6), we used the density functional perturbation theory 
(DFPT) |25| of the linear response of strain type per- 
turbations [26|. The DFPT is based on the systematic 
expansion of the variation expression of the DFT total en- 
ergy in powers of parameters, such as atomic coordinates, 
macroscopic strain and electric field. Tensors x an d c in- 
volve second derivatives of total energy with respect to 
electric field and strain respectively, and tensor e is the 
mixed second derivative of total energy with respect to 
strain and electric field. Here c M „ = with a stress 

and r\ strain. The DFPT calculates C^ v = d fg with E 
the total energy. Under non-zero stress (e.g. pressure), 
Cijki = Cijki — SijOki- The elastic compliance tensor s 
is the reciprocal of tensor c, and the piezoelectric strain 
coefficients di V — Y^u=i e ip. s p.v- We have examined the 
accuracy of c^ v and ei v obtained from the DFPT by com- 
paring with the finite strain method. As summarized in 
Table |U excellent agreement between these two methods 
is achieved. 

We constrained the symmetry to study four phases: 
the paraelectric cubic (C) Pm3m phase, the ferroelec- 
tric tetragonal PAmm, rhombohedral i?3m, and mon- 
oclinic (M) Cm phases. We find that the enthalpy 
(H = E + PV) difference between the T and R phases re- 
duces rapidly with pressure and becomes very small when 
P > 8 GPa. The lowest enthalpy corresponds to the most 
stable phase at that pressure. As displayed in Fig. 1(a) 
and 1(b), there is a T to M phase transition around 9.5 
GPa (V = 56.73 A 3 ), M to R around 11 GPa (V = 56.32 
A 3 ), and R to C around 22 GPa (V = 53.97 A 3 ). The 
M phase for P = 9.5 GPa has its spontaneous polar- 
ization P s along the pseudocubic [uul] direction, where 
u = 0.481. The discovery of a low symmetry M phase is 
completely unexpected in a pure compound like PbTiOa, 
and has only been found in complex solid-solutions like 
PZT [3, PMN-PT 0, etc. 

Because the enthalpy differences among these phases 
are very small, the above calculated pressure-induced 
phase transition sequence needs to be examined by di- 
rectly studying energy derivative properties, such as elas- 
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FIG. 1: Stability of various phases as a function of pressure, 
(a) Enthalpy difference with respect to the C phase for the T 
and R phases, (b) Enthalpy difference with respect to the T 
phase for the R and M phases, (c) Elastic constant C44 and 
en — C12 respectively for the T and R phases, and square of 
the lowest optical phonon frequencies u> 2 for the T [i5(lTO), 
Ai(lTO)], R [-E(ITO), Ai(lTO)], and C [F lu (lTO)] phases, 
(d) Dielectric susceptibility \- 



ticity, phonon frequencies, and diclcctricity. The pressure 
dependence of elastic constants C44 in the T phase and 
en — C12 in the R phase is shown in Fig. 1(c). Note that 
the cubic coordinate system is used for all cases, and elas- 
tic constants Cn = C22 = C33 and C12 = C13 = C23 in the 
R phase. We find that both C44 and cn — C12 approach 
zero around 9 and 11.5 GPa, respectively. Negative C44 
and cn — C12 mean that the T and R phases are unsta- 
ble respectively against shear and tetragonal shear strain, 
indicating phase transitions. The estimated phase tran- 
sition pressures from elasticity anomalies agree well with 
the total energy results. 

In tetragonal PAmm PbTi0 3 , the S(ITO) and 
Ai(lTO) modes originate from the triply degenerate 
Fi u (lTO) mode in the C phase. When the T-to-C phase 
transition occurs at high temperature and ambient pres- 
sure P| or high pressure and room temperature , both 
the iJ(lTO) and j4i(1TO) modes soften simultaneously. 
However, at K we find that pressure induces conden- 
sation only of the E(lTO) mode around P = 10 GPa 
[Fig. 1(c)]. A linear combination of the doubly degen- 
erate i?(lTO) modes gives rise to the transition to the 
M phase. For the R phase, the S(ITO) and Ai(lTO) 
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FIG. 2: Proposed (P, T) phase diagram. It combines theoret- 
ical results (solid circles) at K with experimental data (solid 
squares) at finite temperatures. In this figure we rescaled the 
theoretical pressure for the experimental volume to be zero by 
shifting. The open circles are our guess for these two tricriti- 
cal points, and an ab initio molecular dynamics simulation is 
needed to predict them accurately. 
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modes soften simultaneously at P ~ 22 GPa, going to 
the Fi u (lTO) mode in the C phase. This gives rise to 
the R-to-C phase transition. The linearity of the square 
of mode frequency io 2 vs. pressure close to phase tran- 
sition indicates a Curie- Weiss pressure law. Fig. 1(d) 
shows that the static dielectric susceptibility xn has a 
maximum around 9.5 GPa, indicating a first-order tran- 
sition; while it diverges at P ss 22 GPa, indicating a 
second-order transition, as can be understood from the 
Lyddane-Sachs- Teller (LST) relation. For the T phase 
near transition pressure, xn = X22 c« \/uj 2 [E{\TO)], 
and X33 1/w 2 [^4i(1TO)]. Thus xn increases rapidly 
around 9.5 GPa, whereas X33 does not. The same analy- 
sis can be applied to the R and C phases. The phase tran- 
sition pressures obtained from phonon frequencies and di- 
electric constants also are consistent with the total energy 
results, and it shows that the computed small energy dif- 
ferences between phases are reliable. Our first-principles 
calcalations neglect quantum fluctuations, which could 
reduce the phase transition pressures |27j . 

We have predicted the pressure-induced phase tran- 
sitions of PbTi03 at K. Combining our results with 
experimental finite temperature data, we propose a 
schematic [P, T) phase diagram (Fig. 2), in which the 
intermediate M phase separates the T phase at low pres- 
sures and the R phase at higher pressures. This phase di- 
agram bears remarkable resemblance to the {x, T) phase 
diagram of PZT [J^, where x is the PbTiC-3 composi- 
tion. The narrow M area is the MPB of PbTi0 3 , and it 
serves as a structure bridging the T and R phases in that 
its spontaneous polarization is located in the (110) plane 



FIG. 3: Pressure dependence on piezoelectricity, (a) Sponta- 
neous polarization P s . (b) Piezoelectric stress coefficients eis 
and 633. (c) Piezoelectric strain coefficients dis and 1^33. 



and between the pseudocubic [001] and [111] directions, 
which is very similar to that in PZT. 

Pressure suppresses spontaneous polarization P s [Fig. 
3(a)] and the c-axis strain. Although the magnitude of 
P s changes continuously at the T-M transition pressure, 
Pf is discontinuous, indicating a first-order phase transi- 
tion, consistent with the conclusion drawn from dielectric 
constants. We emphasize that when the T-to-M phase 
transition occurs, P s retains about half of the magnitude 
under zero pressure and c/a — 1.012. These values are 
comparable to those of PMN-PT [l9], PZN-PT [U, and 
PZT 28] in the T phase close to MPB. 

Fig. 3(b) and 3(c) summarize the pressure effect on 
piezoelectricity of PbTiOa. For the T phase, 633 in- 
creases, whereas the piezoelectric strain coefficient 
decreases with pressure. The relatively large e?33 at 
P = —2.2 GPa is due to the LDA overestimation of strain 
at the experimental volume. Both eis and dis increase 
with pressure, and they rise dramatically for P > 4 GPa 
when C44 begins to drop quickly [Fig. 1(c)]. The T phase 
has G?i5 = 615/044, so that d\§ increases even faster than 
eis. The pressure- induced large enhancements of eis and 
d\5 are a result of the sharp reduction of the enthalpy dif- 
ference between the T and R phases, and both eis and 
di^ diverge at the T-to-M phase transition pressure. For 
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the R phase, both eis and di5 are rather small, and 633 
jumps to over 30 C/m 2 , increasing with pressure. d^ has 
huge values near the M phase and remains larger than 
1500 pN/C for a broad range of pressure. The R phase 
under high pressure has (I33 « 633/(011 — C12), and the 
increase of c\\— c\i with pressure [Fig. 1(c)] causes d 33 to 
decrease with pressure, even though 633 rises. Note that 
^33 at 12 GPa is smaller that at 13 GPa, and it could 
be due to numerical uncertainties, since small variations 
of elastic constants C14 and C15 change elastic compli- 
ance s dramatically when en — C12 is small. It is evident 
that the predicted giant piezoelectric effect of PbTi03 
comes from non-collinear polarization rotation since the 
enhancement of di^ in the T phase and d^ in the R 
phase do not occur along the spontaneous polarization 
directions, and the M phase acts as a structural bridge 
between the T and R phases as indicated by pressure 
inducing polarization rotation from [001] to [111] direc- 
tions. 

We have demonstrated that pressure can induce an 
MPB in a simple compound PbTiOs, which is very sim- 
ilar to the composition-induced MPB in complex solid- 
solutions, such as PZT H3], PMN-PT 0, and PZN-PT 
|2fij . It is critical that these ferroelectric systems near 
MPB have distinct phases (different patterns of atomic 
displacements) with very close free energies. Each local 
minimum is shallow and a broad global minimum ex- 
ists. Low symmetry phases can be stabilized between 
two high symmetry phases if both of them become sad- 
dle points, and this is the origin of an MPB. Near an 
MPB, smooth transformation between two high symme- 
try phases with different polarization directions via a low 
symmetry phase results in strong coupling between inter- 
nal degree of freedom and strain. 

In PbTiOs, ferroelectricity arises from the competi- 
tion of short-range repulsions which favor the C phase, 
and Coulomb forces which favor a ferroelectric phase 
0. As pressure is increased, the short-range repul- 
sions promoting paraelectric stability increase faster than 
the Coulomb terms promoting ferroelectric instability. 
Therefore pressure favors the paraelectric phase and sup- 
presses ferroelectricity. Because the tetragonal well- 
depth along [001] reduces faster than the rhombohedral 
wells along [111], under a certain pressure the T and R 
well-depths are equal, and the energy surface connect- 
ing the T and R phase becomes very flat. Higher order 
terms in the energy expression result in an intermediate 
M phase [2^, which has slightly lower energy than both 
of the T and R phases, leading to a T-to-M phase transi- 
tion. Further increasing pressure will make the R phase 
more stable than the M phase. For another important 
ferroelectric, BaTiOs, an M phase is not expected at low 
temperature, because the R well-depth is deeper than the 
T well-depth so that pressure can not induce an M phase 
before both well-depths disappear. 

Pressure-tuning of the structural and piezoelectric 



properties of PbTiOs promises to be an intriguing av- 
enue for experimental study. PbTiOs under pressure 
should be the simplest possible system to study the ba- 
sic physics of the MPB and piezoelectric enhancement 
through polarization rotation. The predicted monoclinic 
Cm phase (without strain and external electric field) and 
the MPB accompanied by giant piezoelectric effects in 
simple perovskite compounds have never yet been re- 
ported. Very recently, Rouquette et al. reported experi- 
mental evidence of polarization rotation induced by pres- 
sure in PZT [3(], El 113 • We hope our predictions will 
stimulate both fundamental and technological interests 
to investigate ferroelectrics under hydrostatic or uniaxial 
pressure, and stress and pressure as means to enhance 
piezoelectricity. 
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